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We present a self-similar model of spherically symmetric collapse of a massless 

O 

scalar field with a parameter p. The black hole formation is explicitly shown to 



occur only in the strong- field implosion of p > 1. The field evolution in the 
critical limit p — > 1 is compared with numerical results found by Choptuik. 



"-j ■ Gravitational collapse is one of interesting phenomena in general relativity. 

In recent years, a large number of studies have been made on the dynamical 
behavior of a massless scalar field, in the context of spherically symmetric gravi- 
tational collapse. It has been shown that for a weak-field implosion scalar waves 
bounce and disperse to infinity, while a strong-field implosion leads to black hole 
formation 0. Furthermore, in numerical calculations, Choptuik found an 
important behavior of a family of solutions which contain the parameter p char- 
acterizing the strength of the scalar field: For black hole formation there exists 
the critical limit p — > p* to give the gravitational mass satisfying a power law 
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Mbh oc| p — p* | 7 with a universal critical exponent 7 ~ 0.37. Such a critical 
behavior may be a general property of gravitational collapse, since imploding 
axisymmetric gravitational waves also show the similar behavior M. 

In this letter our purpose is to present an analytical model of the spherically 
symmetric collapse of scalar waves. We solve the spherically symmetric Einstein 
equations coupled to a massless scalar field ip to find the self-similar solution 
which has a critical parameter p. This model allows us to see explicitly the 
black hole formation and the scalar field dispersion. Our main concern is the 
field evolution in the critical limit p — > 1, which should be compared with the 
numerical results. 

Let us impose a self-similarity on the massless scalar field such that ip(u, v) = 
ip(u/v), where u and v are retarded and advanced times respectively. The line 
element in the spherical symmetric double-null coordinates is given by 

ds 2 = -hdudv + r 2 dQ 2 , (1) 

and we assume the metric to be h(u,v) = h(u/v) and r(u,v) = vf(u/v). For the 
energy-momentum tensor of the massless scalar field defined by 

= ^{d^d^ - i^VM), (2) 

the Einstein equations G^ v = give the unique solution as follows, 

h = 1, (3) 

r 2 = ^{(1 ~p 2 )v 2 - 2vu + u 2 }, (4) 

1 (l-p)v-u 

$ = ±- log , 5 

2 (1 + p)v — u 

where the constant p can be chosen to be positive without loss of generality. In 
this letter we use units such that G = c = 1. 



The space-time is asymptotically flat at the past null infinity u — > — oo, 
where r ~ —u/2. The scalar waves implode from the past null infinity with 
if) ~ ±pv/u ~ ^pv/r. The parameter p turns out to give the amplitude of the 
imploding scalar waves. This solution is useful to analyze the scalar wave col- 
lapse, if a pathological structure due to naked singularities is removed from the 
space-time. Because the Ricci scalar R has the form 

p P 2uv fa\ 
R =^T> ( 6 ) 

the curvature singularity is located at the origin r = 0. Note that the locus of 

r = is given by u = (1 — p)v and u = (1 + p)v in the regions v > and t> < 0, 

respectively. The space-time structure is crucially dependent on the parameter p 

(see Fig. 1). For p > 1 the singularity becomes time-like in the region t> < and 

spacelike in the region v > 0, while for < p < 1 it always becomes timelike. 

(For p — 1 exactly, the curvature singularity at r = in the region v > becomes 

null (i.e., u=0). We do not consider this case, and we will take the critical limit 

p —>■ 1 by keeping the condition p > 1.) The important point is that the apparent 

horizon exists only in the strong-field implosion of p > 1. In fact we obtain the 

expansion 9 for outgoing null rays given by 

e = {{i-p 2 )y-u} 

~ {(l-p)v -u}{(l+p)v + u} 
and the apparent horizon where 9 = is found to be u = (1 — p 2 )v. Thus the 

time-like singularity becomes naked, while the space-like one is hidden by the 

apparent horizon (see Fig. 1). This result will be due to the gravitational mass 

M which can be locally defined in the spherically symmetric system as follows, 

M = ^(l-jrflAr) = -^. (8) 

Near the singularity u = (1 — p)v in the region v > 0, Eq. (|8|) is reduced 
to M ~ p 2 (p - l)v 2 /r. The mass turns out to be negative near the time-like 



singularity (0 < p < 1) and positive near the space-like one (p > 1). Near 
another timelike singularity u — (1 + p)v in the region v < 0, we obtain M < 
irrespective of p. 

In 2-dimensional gravity, the similar relation among the singularity, the local 
gravitational mass and the apparent horizon was pointed out by Hayward in the 
context of the cosmic censorship Note that the 2-dimensional metric can 
have the form ds 2 = r~ 2 dudv with r given by Eq. @ as a special case. This 
may suggest that the 2-dimensional dilation gravity works as a good model of 
the spherically symmetric 4-dimensional gravity. However, in the 4-dimensional 
space-time, the dynamical scattering of imploding scalar waves plays an impor- 
tant role for the evolution of the scalar field. The critical behaviors which will be 
discussed in the following are missed in the 2-dimensional model. 

It is clear that the self-similar solution is not applicable to the region v < 0, 
where the naked timelike singularity u — (1 + p)v appears with M < 0. For a 
realistic evolution we must avoid such a singularity to set some regular initial 
data. Fortunately, this is possible, since the solution requires that there is no 
energy flux across the v = line (i.e. outgoing wave ip u equal to zero.) and the 
mass is equal to zero on it. We can replace the region v < by a flat space-time, 
which is smoothly matched to the region v > (see. Fig 2). For the weak- field 
implosion (0 < p < 1) shown in Fig. 2a, the region corresponding to u > and 
v > is also replaced by a flat metric, and the time-like singularity u — (1 — p)v 
is eliminated from the space-time. (Recall that M = on the u = line.) This 
is the case that the dynamical scattering of the imploding waves can induce a 
complete bounce to the future null infinity v — > oo. For the strong- field implosion 
(p > 1) shown in Fig. 2b, the space-like singularity u — (1 — p)v is hidden by 
apparent horizon. This means a black hole formation by imploding scalar waves. 
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However, our self-similar model requires the imploding amplitude increasing like 
| if) | ~ v and fails to give an asymptotically flat region at the future null infinity 
v — ► oo. (The radius r = p(p 2 — l) l l 2 v/2 of the apparent horizon infinitely 
increases as v — > oo.) Therefore, for gravitational collapse of wave packets shown 
in Fig. 3, the self-similarity must break down in some region denoted by v > vq, 
where outgoing waves dominate incoming ones. Nevertheless the discussions in 
the following suggest that the self-similar solution is approximately valid in a 
strong-field region of < v < v . 

We have found the crucial difference of the space-time structure between the 
strong-field implosion (p > 1) and the weak one (0 < p < 1). Now let us consider 
the critical limit p — > 1 for black hole formation. The numerical calculations of 
wave-packet collapse done by Choptuik claim that the black hole mass approaches 
zero in some critical limit. To discuss such a critical behavior for our self-similar 
model, we must estimate the gravitational mass on the apparent horizon 
instead of the Bondi mass which can be defined only at the future null infinity 
v — > oo. ^From Eq. (|5|) we obtain Mh = p\/p 2 — Iv/A. It is interesting to 
compare with the initial mass Mj = p 2 v/A at the past null infinity u — > — oo. 
The ratio 

M t p ' K } 

may be interpreted as the transmission coefficient of incoming waves propagating 
along a fixed v. For the strong- field limit p >> 1 the amount of dynamical scat- 
tering which generates outgoing waves remains negligibly small before incoming 
waves reach the apparent horizon. For the critical limit p — > 1, however, almost 
all the incoming waves are scattered away, and a small portion of the gener- 
ated outgoing waves can fall into the apparent horizon. (Recall that the surface 
u = (1 — p 2 )v becomes nearly null.) Thus the black hole mass Mh estimated on 
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the apparent horizon vanishes in the limit p — > 1, which gives the power law 

.p — 1 



*v. (10) 

o 

This is slightly different from the power law with the critical exponent 7* ~ 0.37, 
which was numerically evaluated by Choptuik at the future null infinity v — > 00. 
We wish to emphasize that Eq. (|H]) exhibits a mass evolution in some dynamical 
stage of gravitational collapse instead of the asymptotic stage v — > 00. If the self- 
similarity approximately holds in the stage < v < v even for the wave-packet 
collapse shown in Fig. 3, the final mass will be given by Eq. (|i~0"D with vo in place 
of v. Then, to obtain the critical exponent 7* ~ 0.37, the advanced time Vq should 
increase like Vq oc| p — 1 |~ - 13 in the limit p — > 1. For the near-critical evolution 
the self-similar stage may continue for a long time. Though this problem of the 
critical exponent remains unclear within the framework of the present analysis, 
we can point out another critical behavior in relation to the numerical result. 
Now let us introduce the time coordinate t defined by 

l0g(t/r) = " J [ 2-XV P ^ + 4 + P ^ ]dX > (H) 
where X = v/r. The line element can be rewritten into the diagonal form 

ds 2 = -a(X)dt 2 + (3{X)dr 2 + r 2 dQ 2 , (12) 

where the metric coefficients a(X) and /3(X) are 

M h p2X ^-XVp 2 X^ + A + P 2 X 2 i2 ^ 

a{x) = (1 ' v /^ +4 )[ v p^n-p-x ] x 

fl Vp'X- + 4-p 2 X 

eXp(2 i [ 2-X v /^ + 4 + ^ ]rfX) ' (13) 
4 

^ (X) = Vp 2 x 2 + ±Wp 1 x> + 4 - P 2 x) ' ^ 14) 

The scalar field ip also depends on X only, 



^ = ±log-( v /p 2 X 2 + 4+pX). (15) 



We consider the field evolution measured by the coordinates t and r, which are 
related to the variable X through Eq. (|TT|). The variable X satisfies X 2 < X\ = 
4/p 2 (p 2 — 1) outside of the apparent horizon. Then, in the critical limit p — > 1, 
there exists the strong-field region 1 << X 2 << X 2 in a neighborhood of r = 0, 
where Eq. (|TT1) gives the logarithmic relation 



Such a logarithmic evolution was also observed in the numerical calculation of 
wave-packet collapse with a scaling relation. Note that the time "t" propor- 



tional to the centoral proper time which defined by T = f J\ g a (t,0) \dt ~ 



y | p — 1 1 1. For the self-similar solution the scaling relation ?/>(logr — A, logt — 
A) = ip(logr, logt) clearly hold with a continuous parameter A. For the wave- 
packet collapse, however, the scaling parameter A was found to be restricted to 
some discrete values nA* (A* ~ 3.4 and n = 1, 2, 3, . . .). This difference of the 
scaling will not be crucial, if we consider the region where both logr and logt 
are much larger than A*. The apparent horizon can be arbitrarily close to r = 
in the limit p — > 1. Then the relation logr >> A* is satisfied in the strong- field 
region outside the apparent horizon near r = 0, where the spacing of the discrete 
values nA* becomes unimportant to recover approximately the self-similarity. 

In summary, our self-similar model can give a simple picture of black hole 
formation due to imploding scalar waves with the transmission coefficient (Q). 
We obtain the critical behaviors which exhibit the mass evolution fllOD on the 
apparent horizon and the logarithmic evolution (|T6| ) in a neighborhood of r = 0. 
The scaling relation means that the self-similarity may appear even in wave- 
packet collapse with the parameter p sufficiently close to the critical value and 
can represent a dynamical growth of the apparent horizon, though the critical 
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log t — log r. 



(16) 



(p - l)X 2 
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exponent 7* and the scaling constant A* are missed. It is necessary to check the 
self-similar critical behaviors in numerical calculations of wave-packet collapse. 
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Figures 




Fig. 1. 

Penrose diagram for the self-similar solutions covering the whole space-time. 
The scalar waves implode from the past null infinity J~ with the amplitude 
denoted by the parameter p: (a) For the weak-field implosion < p < 1 the 
time-like singularity appears in the negative mass regions of uv > 0. (b) For the 
strong-field implosion p > 1 the apparent horizon (A.H.) is formed to surround 
the space-like singularity in the positive mass region. This space-time has no 
future null infinity J + at v — > oo. 
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Fig. 2. 

Penrose diagram for the self-similar solution, in which the negative mass re- 
gion are replaced by the flat space-time with M — 0. 

(a) For < p < 1 all the imploding waves bounce and disperse to the future 
null infinity J + . (b) For p > 1 we obtain a model of black hole formation due to 
the scalar wave collapse, though the mass infinitely increases as v — > oo. 
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Fig. 3 



Penrose diagram for black hole formation due to wave-packet collapse. The 
region I (v < 0) is the flat space-time, and in the region III (v > v ) outgoing 
waves dominate incoming ones. The self-similarity will remain valid in the region 
II (0 < v < v ). 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/gr-qc/9402017v3 
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